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Abstract: Petrocodon Hance was one of the famous ornamental flowers, but most of them were 
assessed as extremely dangerous (CR), or at least above vulnerable (VU) due to climate 
turbulence and strong interference of human activities. To reconstruct the temporal and spatial 
changes of the potential suitable areas of Petrocodon since the last interglacial period, and explore 
the response of the suitable areas to environmental changes, which provided theoretical guidance 
for the origin of Petrocodon, the study of geographical differentiation, the conservation of unique 
germplasm resources in China, and the development and utilization of gardens. In this paper, 
combined with 120 distribution records and 17 environmental variables, the optimal MaxEnt 
model and geographic information technology (ArcGis) were used to simulate the suitable areas 
and distribution pattern of Petrocodon in China and Indo-China Peninsula. Based on stepwise 
multiple linear regression analysis, redundancy analysis and Monte Carlo test, the dominant 
factors affecting the current geographical distribution of Petrocodon were evaluated. The results 
were as follows: (1) The prediction accuracy of the optimized MaxEnt model was high, and the 
AUC value was greater than 0.96. The potential suitable areas of Petrocodon for the present 
distribution are continuous from southwest China to northern Vietnam, scattered in central and 
southern China and blocky in northern Myanmar, of which the southern Yunnan-Guizhou Plateau 
in China was the most suitable. (2) The environmental variables which restrict the geographical 
distribution of Petrocodon included precipitation of the driest month (bio14), precipitation of the 
warmest quarter (biol8), precipitation of the wettest quarter (biol6), SD of the temperature 
seasonality (bio4), min temperature of the coldest month (bio6) and altitude (alt). (3) Under the 
climate change scenario, the expansion and contraction areas of the suitable habitat of Petrocodon 
were located in the north and northeast of the current potential distribution area, which were 
sensitive areas susceptible to climate change. During the last interglacial period, the suitable area 
of Petrocodon expanded on a large scale, but there was almost no suitable distribution area of 
Petrocodon in the dry and cold environment during the last glacial maximum. After that, with the 
aggravation of climate warming, the suitable habitats of Petrocodon increased rapidly to high 
latitudes, while the suitable habitat in low latitudes decreased. (4) The centroid position of the 
suitable area for Petrocodon migrated northwards from Yongfu County, Guangxi (110.10° E, 
24.69° N) to Chengbu County, Hunan (110.29° E, 26.05° N). To sum up, global warming has a 
positive impact on the potential distribution area of Petrocodon, but extreme warming will cause 
the suitable habitat of Petrocodon to shrink, and the niche of Petrocodon will narrow. Southwest 
China to northern Vietnam, which has the advantage of mature karst landform, may be its main 

refuge. 
Key words: potential suitable areas, climate change, environmental variables, MaxEent model, 
distribution centroid 


全 球 气候 变化 不 断 地 影响 着 物种 的 生活 习性 、 地 理 分 布 、 群 落 组 成 、 生 态 系统 结构 
(Doxford & Freckleton, 2012; 赵 泽 芳 等 ，2016; Matias et al., 2017) ， 导 致 某 些 对 生境 专属 
性 要 求 较 强 的 植物 迁移 无 路 ， 加 速 物种 灭绝 ， 同 时 也 改变 物候 节律 引发 生态 紊乱 (Sandel et 
al., 2011; Li & Chen, 2014) 。 气 候 变 化 引起 地 表 温 度 和 降水 的 重新 分 配 ， 而 地 表 温 度 和 降水 
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决定 了 植被 的 形成 和 演化 (Lambert et al., 2010; Wang et al., 2014) ， 进 而 显著 影响 扩散 能 
E DREHEN A SE E EF j C Walther et al., 2002; Hewitt, 2004; Svenning & Skov, 
2007) 。 能 量 是 植物 生命 活动 的 基础 和 能 量 来 源 ， 水 分 是 植物 生理 活动 的 来 源 和 基本 成 分 
(Fang, 1991) 。 因 此 ， 关 注 过 去 和 未 来 物种 对 气候 变化 的 响应 ， 不 仅 有 助 于 了 解 物种 形成 
的 历史 原因 和 地 理 分 布 的 变化 ， 而 且 有 助 于 制定 科学 的 种 质 资源 管理 策略 。 
根据 当前 文献 报道 ， 石 山 莒 苔 属 CPetrocodon Hance) 主要 分 布 于 中 国 及 中 南 半 岛 的 喀 
斯 特地 貌 区 域 ( 卢 永 彬 等 , 2017) , 共计 47 种 ( 含 1 变种), 其 中 中 国 报道 44 种 (GRC, 2022; 
温 放 等 , 2022b) 。 该 属 植物 花 部 变异 几乎 涵盖 了 所 有 的 昔 曹 苦 科 的 花 部 特征 〈 葛 玉 珍 , 2012; 
卢 永 彬 等 ，2017) ， 体 态 多 样 ， 花 色 艳 丽 ， 花 冠 局 部 还 常 兼 有 不 同 颜色 的 色 晕 、 和 条纹、 斑点 
或 网 纹 ， 具 有 良好 的 产业 开发 应 用 前 景 〈 卢 永 梢 等 ，2017) 。 该 属 植物 分 布 于 阴 湿 的 喀斯特 
地 貌 (少量 种 类 见于 丹 霞 地 貌 ) 的 峡谷 、 石 缝 、 岩 壁 或 溶洞 口 弱 光 带 ， 对 石灰 岩 钙化 土 的 基 
质 具 有 专 一 性 ， 局 限 分 布 于 一 个 大 区 域内 生态 环境 最 为 优越 的 位 置 〈 温 放 ，2021) 。 并 且 ， 
该 属 植物 多 数 物种 数量 少 ,分布 范围 窄 ， 对 环境 的 要 求 严格 并 且 环 境 适 应 脆弱 ， 是 评价 当地 
生态 环境 优 劣 的 重要 植物 类 群 之 一 〈 韦 角 刚 ，2018; 辛 子 兵 等 ，2019) 。 一旦 它们 赖 以 生存 
的 生态 位 发 生变 化 ， 就 意味 着 这 个 物种 或 这 个 地 区 某 个 群体 的 消失 《〈 葛 玉 珍 ，2012; 温 放 ， 
2021) 。 值 得 注意 的 是 ， 由 于 气候 动荡 和 人 类 活动 的 强烈 干扰 ， 近 年 发 表 的 类 群 中 除 小 黄花 
^ai EZ luteoflorus Lei Cai & F.Wen) (Fan et al., 2020 ) 外 ， 其 余 至 少 也 是 易 危 CVU) 
以 上 ， 多 数 被 评估 为 极 危 CCR) (Li etal., 2019; Su et al., 2019; Zhang et al., 2019; Fan et al., 
2020; Li et al., 2020; Xin et al., 2021; Nong et al., 2021) 。 
© 目前 , WHE ARNE HEARN RR, ARK A MBE IK BEYE 
= ffi CBB, 2012; 温 放 等 ，2022a) . WALI SRI CEA ^E RIE RA EA FIG 
Pap 未 进行 系统 和 深入 的 研究 ,诸多 价值 还 未 挖掘 出 来 ， 该 属 植物 就 已 经 处 于 极 危 的 状态 。 虽 已 
开展 大 量 野外 调研 和 濒危 现状 的 再 评估 , 但 野外 调查 获得 的 数据 和 濒危 等 级 不 足以 充分 反映 
该 属 植物 的 整体 分 布 格局 。Liu 等 (2017) 认为 石山 芒 蔡 属 物种 分 布 与 地 理 坏 境 有 着 直接 关 
系 ,受到 诸多 环境 变量 的 制约 ,石山 芒 蔡 属 是 在 特定 地 区 或 某 种 局 部 特殊 生境 下 形成 的 物种 
具有 空间 信息 特征 ( 温 放 等 ，2022b) 。 因 此 ， 重 建 不 同 气候 条 件 下 石山 苔 苔 属 在 中 国 及 中 
南 半 岛 的 分 布 格局 ， 分 析 其 生态 适宜 性 意义 重大 。 
= 物种 能 够 在 特定 地 区 或 某 种 局 部 特殊 生境 中 生存 和 繁衍 后 代 的 所 有 条 件 总 和 称 为 生态 
一 位 TAPES, 2013) 。 生 态 位 的 研究 可 以 通过 生态 模拟 实现 〈 朱 耿 乎 等 , 2013; Ahmed et al., 
€) 20150 。 生 态 模拟 是 通过 收集 物种 已 知 分 布 数据 和 生物 环境 变量 ， 再 结合 统计 学 的 运算 方法 
算 研 究 物 种 的 生态 需求 , 并 将 运算 结果 投射 到 不 同 的 时 空中 用 以 预测 物种 的 潜在 地 理 分 布 
( 朱 耿 平等 ，2013) > A, MARAMA! CMaxEnO 模拟 植物 潜在 地 理 分 布 已 经 成 为 
1 研究 趋势 ， 并 得 到 较 好 的 研究 结果 (Higgins etal., 2012; Cory et al., 2013; Higginset al., 
2020; Zhao et al., 2021) 。 特 别 是 将 MaxEnt. ArcGis. R 语言 及 ENMTools Z (44044 E 
其 他 物种 分 布 模型 比较 时 ， 结 果 更 接近 于 已 实现 的 生态 位 (Higgins etal., 2012) ， 并 且 预 测 
的 结果 便于 分 析 CAhmed et al., 2015) , 准确 度 也 相对 较 高 (Phillips & Dudik, 2008; Elith et al., 
2011) ， 尤 其 在 研究 物种 分 布 广泛 和 位 置信 息 不 确定 中 有 很 大 的 优势 CFarashi et al., 2013) 。 
因此 , MaxEnt 模型 是 目前 使 用 最 为 广泛 的 生态 位 模型 (Ahmed et al., 2015; Barbosa & Schneck, 
2015; Vaz et al., 2015) > MaxEnt 模型 有 助 于 推断 物种 已 知 分 布 范 围 内 的 潜在 分 布 ， 但 由 于 
生物 相互 作用 、 扩 散 限制 和 在 不 适 环境 中 的 持久 性 〈Higgins etal., 2020) ， 可 能 会 限制 模型 
的 外 推 能 力 ， 对 物种 已 知 分 布 范 围 外 的 分 布 预 测 存 在 局 限 性 (Sillero, 2011) 。 因 此 ， 本 万 
究 选用 MaxEnt 模型 来 模拟 不 同 气候 条 件 下 石山 营 谷 属 在 中 国 及 中 南 半 岛 的 潜在 分 布 区 。 
当前 , 关于 石山 营 昔 属 适 生 区 及 其 与 气候 变化 关系 的 研究 还 很 薄弱 , 该 研究 资料 的 缺乏 
不 利于 该 属 植物 的 价值 评估 和 资源 管理 。 本 研究 以 中 国 和 中 南 半 岛 为 研究 区 域 ， 选 取 17 个 
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环境 变量 ， 依 托 石 山 曹 营 属 地 理 分 布 信息 ， 采 用 生态 学 模型 和 地 理 信息 技术 的 方法 ， 利 用 
MaxEnt 模型 、ArcGis、R 语言 及 ENMTools 多 软件 相 结 合 ， 通 过 模拟 预测 石山 昔 苔 属 的 潜 
在 分 布 区 和 适宜 生境 , 拟 探讨 以 下 问题 ，(1) 预测 石山 营 苔 属 适 生 区 的 时 空 分 布 格局 ， (2) 
探索 限制 当前 石山 营 苔 属地 理 分 布 的 重要 环境 变量 , 综合 分 析 该 属 植 物 地 理 分 布 与 环境 变量 
之 间 的 关系 ; G) 提出 保护 和 开发 利用 石山 营 苔 属 植物 的 科学 依据 。 
1 材料 和 方法 
1.1 研究 区 概况 
石山 营 苔 属 集中 分 布 中国 及 中 南 半 岛 ， 分 布 范围 为 98.8° 一 117.2° E、10.7° 一 33.3° N; 
分 布 海拔 为 120—1 700 mm。 从 地 形 分 布 来 看 ， 研 究 区 呈现 北 高 南 低 的 格局 〈 钱 满 亿 ，2021 )， 
是 亚洲 热带 的 北 缘 ， 也 是 南亚 季风 、 东 亚 季 风 和 西北 太平 洋 季风 的 受 加 影响 区 ， 正 是 高 山大 
河 与 喀斯特 地 貌 集 中 分 布 的 中 国 西南 及 中 南 半 岛 一 带 ( 姜 超 等 ， 2017) 。 在 当地 复杂 地 形 地 
貌 及 多 种 气候 条 件 的 影响 下 ,形成 了 多 种 多 样 的 局 域 性 小 气候 和 地 带 性 植被 (Zhu, 2013; Ren, 
2015) ， 是 热带 植物 适应 分 化 形成 大 量 区 域 特有 物种 的 一 个 热点 [X (Myers et al., 2000) . 
1.2 数据 来 源 
1.2.1 物种 分 布 的 收集 与 整理 
石山 曹 苔 属 的 分 布 记录 数据 主要 来 源 : CL 课题 组 野外 调查 ; (D 植物 标本 和 数据 库 。 
包括 :中 国 数字 植物 标本 馆 CVH (https://www.cvh.ac.cn/) 、 全 球 生 物 多 样 性 信息 网 络 GBIF 
Chttps://www.gbif.org/) 、 世 界 植物 在 线 POWO Chttps;//powo.science.kew.org/) 。 (3) X 
献 记载 。 为 了 避免 因 记 录 点 之 间 的 距离 过 近 引 起 过 度 拟 合 。 首 先 , 将 收集 的 分 布 数 据 按照 物 
Phá, ARE QO , RE (YO , JEU C csv) 文件 类 型 进行 保存 ， 共 计 得 到 180 个 分 布 数 
据 。 其 次 ， 为 了 保证 每 一 个 物种 有 1 分 布点 ， 将 分 布 记 录 小 于 3 的 物种 直接 用 于 模型 预测 ， 
再 将 分 布 记录 大 于 3 的 种 进入 下 一 步 筛选 。 紧 接着 ,以 10 km 范围 内 的 多 个 样本 点 只 保留 一 
个 点 的 原则 ， 在 ENMTools 软件 中 利用 除 元 余数 据 (trim duplicate occurrences) 工具 删除 10 
km X 10 km 栅 格 内 重复 记录 数据 (Bian & Shi, 2019) . ALA Limi, RAIK 120 份 
石山 营 苔 属 样本 分 布 信息 。 
WW 1.2.2 环境 变量 来 源 及 筛选 
一 该 文 研究 共计 35 个 环境 变量 ( 表 1 ) 。 其 中 ，19 个 生物 气候 变量 图 层 来 自 世 界 气象 
Oo 数据 库 Chttps://www.worldclim.org/) ， 气 候 系统 模型 为 CCSM4， 该 气候 模型 应 用 于 气候 模 
拟 具 有 较 大 优势 〈 肖 建华 等 ，2021) 。 气 候 数据 共计 6 个 时 期 ， 分 别 为 末次 间 冰 期 [last 
interglacial (LIG) , 距 今 约 120 000 Æ 140 000 年 ]、 末 次 冰 盛 期 [last glacial maximum (LGM), 
KA) 22 000 年 前 ]、 全 新 世 中 期 [mid-holocene (MH) ， 大 约 6 000 年 前 ]、 当 前 时 期 (1950 
一 2000 年 ) 2050 Æ (2041—2060 年 ) 和 2070 4E. (2061—2080 年 ) 。13 个 土壤 变量 数据 
来 JR 农 组 织 ( 联合国 粮食 及 农业 组 织 ) 
( https://www.fao.org/land-water/databases-and-software/hwsd/en/ ) 的 世界 土壤 数据 库 
( HWSDO 。 地 形变 量 中 的 海拔 变量 来 世界 气象 数据 库 
Chttps://www.worldclim.org/data/worldclim21.html) ， 坡 度 变量 和 坡 向 变量 通过 海拔 变量 在 
ArcGis 中 利用 表面 分 析 工 具 获 得 ,矢量 地 图 数据 来 自 自 然 资源 部 (http://bzdt.ch.mnr.gov.cn/)。 
在 实验 中 , 为 了 预测 模型 的 顺利 进行 ， 需 要 将 所 有 环境 变量 利用 统一 行政 边界 进行 裁剪 ， 并 
有 具 有 统一 的 栅 格 数据 像 元 大 小 (0.008 3, 0.0083) 、 地 理 坐 标 系 CWGS 1984) 以 及 投影 坐 
标 系 (WGS 1984 EASE Grid Global) 。 
在 进行 物种 地 理 分 布 模拟 时 , 为 了 避免 各 环境 变量 多 重 共 线 性 而 导致 模型 预测 结果 的 过 
度 拟 合 问 题 ， 需 要 预先 运行 MaxEnt 模型 ， 通 过 Jackknife JINA) 检验 环境 变量 的 重要 性 
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及 贡献 率 , 剔除 贡献 率 为 


中 ， 利 用 提取 工具 中 的 sample 工具 进行 采样 ， 把 点 
中 分 析 石 山 划 车 属 分 布 与 环境 变量 之 | 
再 利用 ENMTools 对 环境 变量 进行 上 尔 森 相关 愧 
时 , 删 减 对 预测 贡献 率 较 小 的 


1 ), 当 相 关系 数 大 于 0.8 


0 的 变量 (Zhu & Qiao, 2016) 。 再 将 分 布 数据 和 环境 变量 导入 ArcGis 


带 图 层 提 取 为 excel 表格 ， 然 后 在 spss 
司 的 相关 性 ， 相 关系 数 大 于 0.95 的 环境 变量 得 以 保留 。 

EZ NT, 并 在 R 中 绘制 皮尔 森 相 关 性 热 图 (图 
变量 。 经 筛选 得 到 9 个 气候 变量 (bio2、 


bio3 bio4. bio6. biol4, biol5, biol6. biol8. biol9) 、5 个 土壤 变量 (t texture、t_caco3、 
和 3 个 地 形 地 貌 变量 Calt, aspect. slope) 。 
dé 1 所 有 环境 变量 的 描述 


Table 1 Description of all environmental variables 


t oc. t sand, s gravel) 


类 型 Type fi Description 
biol 年 平均 大 气温 度 Mean annual temperature 
ae Mean diurnal range [Mean of monthly( max 
bio2 昼夜 温差 月 均值 
temperature-min temperature) 
bio3 昼夜 温差 与 年 温差 比值 Isothermality ( bio2 / bio7 X 100) 
bio4 温度 变化 方差 SD of temperature seasonality 
bio5 最 热 月 份 最 高 温 Max temperature of the warmest month 
bio6 最 冷 月 份 最 低温 Min temperature of the coldest month 
bio7 Em FEAR P Temperature annual range ( bio5 — bio6) 
bio8 最 湿 季 度 平均 温度 Mean temperature of the wettest quarter 
气候 Climate bio9 最 干 季度 平均 温度 Mean temperature of the driest quarter 
biol0 最 热 季 度 平均 温度 Mean temperature of the warmest quarter 
bioll 最 冷 季 度 平均 温度 Mean temperature of the coldest quarter 
bio12 年 均 降雨 量 Annual precipitation 
biol3 最 湿 月 降雨 量 Precipitation of the wettest month 
biol4 最 干 月 降雨 量 Precipitation of the driest month 
biol5 湿度 变化 方差 SD of humidity seasonality 
biol6 最 湿 季 降雨 量 Precipitation of the wettest quarter 
biol7 最 干 季 降雨 量 Precipitation of the driest quarter 
biol8 度 Mean precipitation of the warmest quarter 
biol9 Mean precipitation of the coldest quarter 
t_ph Topsoil pH (H20) 
t_clay 顶层 土壤 粘土 含量 Topsoil clay fraction 
t_sand 页 层 土 壤 沙 含量 Topsoil sand fraction 
t_texture 顶层 土壤 质地 Topsoil calcium carbonate 
t caco3 顶层 土壤 碳酸 盐 含量 Topsoil calcium carbonate 
t_gravel 顶层 土壤 碎 石 体积 百分比 Topsoil gravel content 
土壤 Soil t_oc 顶层 土壤 有 机 碳 含量 Topsoil organic carbon 
s ph 底层 土壤 酸碱度 Subsoil pH (H20) 
s clay 底层 土壤 粘土 含量 Subsoil clay fraction 
s sand 底层 土壤 沙 含量 Subsoil sand fraction 
S caco3 底层 土壤 碳酸 盐 含量 Subsoil calcium carbonate 
s gravel 底层 土壤 碎 石 体积 百分比 Subsoil gravel content 
S oc 底层 土壤 有 机 碳 含量 Subsoil organic carbon 


地 形 地 貌 Land alt 


Altitude 
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类 型 Type 描述 Description 
slope 坡度 Slope 
aspect 坡 向 Aspect 
2 o p o p 
2 6 8 >g E 2 
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bio4 e 00000000000 [EOD m 
bio5 000000 0000 o o : 
bio6 © 00000000000 o 
bio7 © 0000000000 Goce 
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S OC Q 
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环境 变量 缩写 详 见 表 1. 
See table 1 for abbreviations of environmental variables. 
1 MAARE YB RD BE 


Fig. 1 Pearson correlation matrix of all environmental variables 


1.3 模型 参数 优化 

MaxEnt 模型 参数 设置 的 复杂 程度 影响 模型 的 转移 能 力 Merrow et al., 2013; Qiao et al., 
2015) ， 而 转移 能 力 低 的 模型 会 造成 预测 结果 不 可 靠 或 很 难 解释 〈 朱 耿 平 和 乔 慧 捷 ，2016 )。 
研究 表明 ，MaxEnt 模型 的 复杂 度 与 正则 化 乘 数 (regularization multiplier, RM) 和 特征 组 合 

(feature combination, FC) 参数 密切 相关 (Muscarella et al., 2015; 朱 耿 平和 乔 慧 捷 ，2016 ) 。 
可 以 通过 调用 Kuenm 数据 包 调控 RM 和 FC, 分 析 各 模型 参数 的 复杂 程度 , 选取 复杂 度 低 而 
转移 能 力 强 的 模型 参数 〈 朱 耿 平 和 乔 慧 捷 ，2016) . MaxEnt 模型 默认 参数 有 5 种 特征 ， 分 
别 为 线性 〈linear-L) 、 二 次 型 〈quadratic-Q ) 、 片 段 化 Chinge-H) 、 乘 积 型 (product-P) 
MBETE Cthreshold-T) ; RM 的 值 为 1。 

在 R 中 ， 调 用 Kuenm 数据 包 优化 MaxEnt 模型 (Cobos et al., 2019) ， 将 正则 化 乘 数 设 
AY 0.1—4, 每 次 间隔 0.1, 一 共 40 种 调控 倍 频 ; 采用 6 个 特征 组 合 , 即 : L、LQ、H、LQH、 
LQHP、LQHPT RICE ATE RUE, 2016) > Kuenm 数据 包 对 上 述 31*40 种 参数 组 合 进行 测 
试 。 当 满足 pval pROC 显著 ， 遗 漏 率 (omission rates) <5%, W AICc-2, delta AICc-0, 


则 特征 组 合 模型 为 最 优 (Phillips et al., 2017; Cobos et al., 2019; Zhuo et al., 2020) 。 


1.4 模型 模拟 


将 经 过 筛选 得 到 的 17 个 环境 变革 
数 设置 如 下 : 随机 测试 集 为 25%; 
变量 的 重要 性 和 贡献 率 ; fi 


复 运 行 10 次 ,最终 输出 结果 文件 为 10 次 的 平均 值 。 除 此 之 外 ， 其 他 均 保 持 默认 参数 。 利 用 
ArcGis 10.4 对 模型 输出 的 栅 格 地 图 进行 再 次 编辑 ， 
从 0 到 1。 然 后 ， 用 自然 分 割 法 (Jenks〉 对 分 布 预测 图 进行 重 分 类 ， 将 潜在 分 布 区 


个 等 级 ， 再 利用 手动 法 输入 


E 


2 


[X (0.3~0.5) 、 最 适 生 区 (05~1) 。 在 实验 中 ， 高 适 生 区 和 最 适 生 


1.5 石山 鞋 苦 属 自然 分 布 对 环境 变量 的 响应 关系 


F (RDA) ， 并 采 


j 刀 切 法 评估 环 | 


EI 120 个 分 布 记录 分 别 导 入 MaxEnt 3.4.1 中 。 模 型 参 
由 环境 响应 曲线 和 预测 分 布 图 ， 利 
出 格式 为 Logistic; RM 和 FC 根据 模型 优化 结果 设置 。 模 型 习 


am 


pd 


每 个 栅 格 像 元 的 值 表示 存在 概率 ， 范 转 


WE: MEEK 〈0 一 0.1) 、 低 适 生 区 〈0.1 一 0.3) 、 高 适 入 
广 定 义 为 总 适 生 区 。 


分 为 4 


利用 SPSS 20.0 软件 进行 描述 性 统计 来 分 析 石 山 童 苦 属 适 生 区 的 气候 特征 ， 并 采用 变 
RAB (CV) 来 衡量 各 气候 变量 的 变异 程度 (Fang & Yoda, 1991) 。 最 后 ， 分 别 将 17 个 环 
Bak ERI E ESTE EBERT 


Hj UE aet AARET JURA 


回归 分 析 ， 并 在 CANOCO 5.0 软件 中 采 
蒙特 卡 洛 (Monte-Carlo) 检验 量化 各 


气候 变量 对 石山 曹 苦 属 地 理 分 布 的 贡献 率 , 以 此 来 判断 影响 石山 营 苔 属 自 然 分 布 的 主导 气候 


变量 。 
1.6 适宜 分 布 区 空间 格局 


物种 存在 概率 值 过 0.30 的 空间 单元 定义 为 石山 营 苔 属 适宜 区 域 ， 物 种 存在 概率 值 <0.30 


的 空间 单元 被 认为 是 不 适宜 


分 布 存在 /不 存在 (0，1) HB 
示 不 存在 的 值 (0〉。 以 矩阵 表 0、1 为 基 而 
属 适 宜 分 布 区 空间 格局 的 变化 。 定 义 了 4 类 适宜 区 变化 : 新 增 适 宜 区 、 缩 减 适宜 区 、 保 留 适 
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x 域 。 建立 了 当前 、 过 去 和 未 来 气候 情景 下 石山 营 苔 属 潜在 地 理 
ERE 1) ， 不 适宜 生境 被 赋予 表 


上 ， 进 一 步 分 析 了 过 去 和 未 来 气候 情景 下 石 


al 
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1.7 质心 转移 分 析 


首先 , 利用 质心 位 置 的 变化 来 反映 物 


而 得 。 


WE 


:适宜 生境 的 迁移 方向 和 距离 (Zhao et al., 2021) 。 


其 次 ， 利 用 及 软件 中 的 Kuenm 包 计 算 了 末次 间 冰 期 、 末 次 冰 盛 期 、 全 新 世 中 期 、 当 前 和 未 


来 适宜 区 域 的 空间 变化 和 适宜 区 域 的 几 
2018) 。 最 后 ， 我 们 使 用 SDMs 1 
心 位 置 ， 计 算 质 心经 纬度 坐标 并 进 


2 结果 与 分 析 
2.1 模型 优化 及 准确 性 评价 


基于 石山 曹 苦 属 的 120 个 分 布点 和 17 个 环境 变量 ， 利 


[有 具 跟踪 质心 的 经 纬度 ， 以 检查 不 同时 期 石 ] 


步 得 出 迁移 昌 


BRS (Cong etal., 2020) 。 


EE 


可 质心 位 置 (Bateman et al., 2016; Laurent et al., 


属 的 质 


] MaxEnt 模型 模拟 和 预测 该 属 


植物 潜在 的 适宜 分 布 区 。 根 据 模型 优化 结果 ( 表 2 )， 满足 pval_pROC=0 具有 统计 显著 性 、 
遗漏 率 =0.0333% 小 于 5%, 也 小 于 MaxEnt 模 型 默认 值 . W_AICc=1 小 于 2, 与 其 对 应 的 RM=0.1 


fll FC=L. 
d 2 Kuenm 包 生 成 的 MaxEnt 模型 的 评价 指标 
Table 2 Evaluation metrics of MaxEnt model generated by Kuenm 
类 型 特征 组 合 正则 化 乘 数 遗漏 率 
Tos Feature combination Regularization multiplier etur 
(FC) (RM) 
默认 值 Default value LQH 1 0.192 3 
优化 值 Optimized value L 0.1 0.033 3 


受 试 者 工作 特性 曲线 (receiver operating characteristic curve, ROC) 面积 在 0.5~1.0. 
一 般 认 为 AUC 值 在 0.5~0.7 时 ,模型 预测 可 靠 性 较 低 ; 在 AUC 0.7—0.9 时 ,模型 预测 可 靠 
性 中 等 ，AUC 大 于 0.9 时 ， 模 型 预测 可 靠 性 最 高 (Hosmer & Lemeshow, 2000) 。 我 们 的 古 
究 结 果 表 现 为 各 时 期 AUC 值 均 >0.9 (图 2 ) ， 表 明 模 型 预测 可 信 度 高 。 
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49-5 E Specifity 特异 度 Specifity 特异 度 Specifity 
A. 末次 间 冰 期 (LIG， 距 今 约 120 000—140 000 年 ) ; B. 末次 冰 盛 期 LGM， 大约 22 000 年 前 ) ; C. 全 
折 世 中 期 (MH， 大 约 6 000 年 前 ) ; D. 当前 (1950—2000 年 ) ;，E. 2050 年 (2041 一 2060 4E) ; F.2070 
年 (2061 一 2080 年 ) 。 
A. Last interglacial (LIG, about 120 000-140 000 years BP); B. Last glacial maximum (LGM, about 22 000 years 
ago); C. Mid-holocene (MH, about 6 000 years ago); D. Current (1950—2000); E. 2050 (average for 2041—2060); 
F. 2070 (average for 2061—2080). 


> 


2 MaxEnt 模型 的 ROC 检验 
Fig. 2 ROC test of MaxEnt model 
2.2 限制 当前 潜在 分 布 区 格局 的 环境 变量 
石山 营 苔 属 存在 概率 与 环境 变量 之 间 的 逐步 线性 回归 关系 如 表 3 所 示 : eG EE 
属 存 在 概率 的 水 分 变量 有 最 干 月 降雨 量 (bio14) 、 最 冷 季 度 平均 降雨 量 (bio19) 、 最 热 
度 平均 降雨 量 (bio18) 、 最 湿 季 降雨 量 (bio16) ; 能量 变量 有 昼夜 温差 月 均值 (bio2) 、 
昼夜 温差 与 年 温差 比值 (bio3 ) 、 温 度 变 化 方差 (bio4) 、 最 热 月份 最 高 温 (bio5) ; HUE 
变量 有 海拔 (alt) 、 坡 向 (aspect) ; 土壤 变量 有 表层 土壤 碳酸 钙 (t caco3) 、 底 层 土壤 砾 
含量 (s gravel) 、 表 层 土壤 砂 含量 (t sand) 、 表 层 土 壤 质 地 (t texture) 。 以 上 主导 变 
量 解释 率 高 达 98.1%， 且 均 达 到 极 显 著 水 平 (P<0.01)。 坡 度 (slope) 、 表 层 有 机 碳 (t_oc)、 
湿度 变化 方差 (bio15) 3 个 气候 变量 在 拟 合 过 程 中 被 剔除 ， 表明 这 3 SUR LE E 
属 的 存在 概率 影响 较 小 。 
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表 3 石山 童 苦 属 分 布 区 环境 变量 的 描述 性 统计 


Table 3 Descriptive statistics of environmental variables in the distribution regions of 


Petrocodon 
环境 变量 相关 系数 Fi 显著 性 


Significant differences 


(S YEUX AR 


=> 


Environmental variable Correlation coefficient F value 


(RA (PY 


Durbin-Watson 


环境 变量 相关 系数 pik 显著 性 德 宾 沃 森 
E Correlation coefficient sige Significant differences ERE 
R (P) 

biol4 0.981 614.805 <0.01 
alt 0.98 615.079 <0.01 
aspect 0.973 506.831 <0.01 
biol9 0.967 444.457 <0.01 
bio3 0.96 406.324 <0.01 
t_caco3 0.953 382.767 <0.01 

s_gravel 0.945 366.896 <0.01 gi 
bio4 0.934 346.409 «0.01 
t sand 0.916 313.274 «0.01 
biol6 0.897 301.24 «0.01 
bio6 0.87 289.865 «0.01 
biol8 0.824 273.265 <0.01 
bio2 0.766 287.581 <0.01 
t_texture 0.632 303.394 <0.01 


为 了 进 


后 ， 利 


TE: 环境 变量 的 含义 见 表 1。 下 同 。 


Note: See Table 1 for meaning of environmental variable. The same below. 


A PRE A LL E EE Js 2 [n] AE ES PINE E 
量 和 120 DREHER, TE CANOCO 5.0 软件 
蒙特 卡 洛 检验 14 个 环 : 


WiZ EB AE gc Ase UA LLL ER 2S Je 3.3 
的 有 6 个 环境 变量 ， 分 别 为 最 干 月 
(biol6) 、 温 度 变 化 方差 (bio4) 、 最 冷 月 份 最 低温 (bio6) 、 


`Œ (P<0.005) 


(biol8) 、 最 湿 季 降雨 量 


海拔 Cal 。 其 中 ， 


解释 率 和 贡献 率 分 别 为 10.7% 和 11.7%; E 


降雨 量 


的 关系 ， 将 以 上 筛选 得 到 14 个 环境 变 
进行 见 余 分 书 
漠 变 量 对 地 理 分 布 的 解释 力 排序 
分 布 格局 的 主导 环境 变量 。 其 


F (CRDA) ， 经 前 向 选择 
#4). 4 P<0.005 时 ， 
， 贡 献 率 达到 极 显 著 水 
& (biol4) 、 最 热 季 度 平 均 降 雨量 


能 量变 量 bio4 和 bio6 累计 解释 率 和 累计 贡献 率 最 高 ， 分 别 为 75.5%、 
79.4%; 水 分 变量 bio14、bio18、bio16 累计 解释 率 和 累计 贡献 率 为 7.4% 和 8%; 地 形变 量 的 


寞 变量 的 PP 值 均 大 于 0.05， 未 达到 显著 水 平 。 因 


此 ， 影 响 石 山 苔 苔 属地 理 分 布 的 主导 环境 变量 有 bio4、alt、bio14、bio16、bio18、bio6。 各 


类 环境 变量 对 石山 童 苦 属 地 理 分 布 格局 的 影 


自力 排序 为 能 量变 量 > 地 形变 量 > 水 分 变量 。 


表 4 环境 变量 的 蒙特 卡 洛 检验 结果 
Table 4 Monte Carlo test results of environmental variables 
环境 变量 解释 率 贡献 率 Fit 显著 性 
ne ae ee Explains rate Contribution rate an Significant differences 
(%) (%) (P) 
bio4 71.1 77.9 436 0.002 
alt 10.7 11.7 103 0.002 
biol4 3.9 4.2 47.3 0.002 
biol6 1.1 1.2 14.7 0.002 
bio18 2.4 2.6 38.5 0.002 
bio6 4.4 1.5 30.5 0.002 
bio2 0.3 0.3 4.7 0.02 
aspect 0.4 0.4 6.7 0.01 
slope 0.2 0.2 3.9 0.062 


环境 变量 解释 率 贡献 率 Ft 显著 性 
uL UM Explains rate Contribution rate uis Significant differences 
(%) (70) (P) 

bio3 0.2 0.2 3.4 0.072 
biol9 0.1 0.1 2.4 0.104 
t_caco3 0.1 0.1 23 0.122 
t texture 0.2 0.2 2.8 0.092 
biol5 0.1 0.1 24 0.134 
t oc «0.1 «0.1 0.8 0.412 
s gravel «0.1 «0.1 0.3 0.674 
t sand «0.1 «0.1 0.2 0.778 


对 6 个 显著 影响 石山 营 苔 


3) ， 结 果 表 明 前 2 个 3 


分 轴 的 


属地 理 分 布 的 气候 变量 与 地 到 


解释 率 分 别 为 73.32% 和 90.2796. 


分 布 数据 进行 RDA 排序 (图 


其 中 ， 经 度 (X) 5j biol4. 


bio4 呈正 相关 关系 ， 与 alt、bio6、bio18 呈 负 相关 关系 ， 与 bio16 不 相关 ; AL (YO 与 bio4、 
biol4 呈正 相关 关系 ， 与 bio16、bio6 呈 负 相关 关系 ， 与 alt、bio18 不 相关 。 
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X. 经 度 ; Y. AE, + RRO WES RHE 
和 纬度 ， 空 心 箭 头 表 示 自 变量 的 气候 变量 ; 


直角 表示 不 相关 性 ) 。 


夹 角 大 小 表示 相关 忆 


轴 1 Axis 1 (90.27%) 


E 《锐角 表示 正 相 关 性 ; 钝 角 表 示 负 相关 性 ; 


1.0 


X. Longitude; Y. Latitude; * indicates locations of Petrocodon. See table 1 for meaning of climatic variables. The 


solid arrows indicate the longitude and latitude of the dependent variables, the hollow arrows indicate the climate 


variables of the independent variables; The angle indicates the correlation (acute angle indicates positive 


correlation; obtuse angle indicates negative correlation; right angle indicates irrelevance). 


eh Se 


到 3 石山 


ni 


EA EB 


的 地 理 分 布 与 环境 变量 的 RDA 排序 


Fig. 3 RDA ordination of geographical distribution and environmental variables of Petrocodon 


为 了 解释 6 个 主导 环境 变量 与 石山 


释 物 种 存在 概率 和 环境 变量 之 i 


司 的 关系 (EAD. 


属地 理 分 布 的 关系 , 通过 绘制 环境 响应 曲线 来 解 


其 中 ， 适 


宜 区 的 温度 变化 方差 Cbio4) 


Im 


的 范围 为 98 一 920、 最 热 季 度 平 均 降 雨量 (bio18) 的 范围 是 185 一 1 207 mm. HFK (alt) 的 
范围 是 70—2 084 m、 最 干 月 降雨 量 (bio14) 的 范围 是 1 一 54 mm、 最 湿 季 降雨 量 (bio16) 
的 范围 是 406—1207 mm、 最 冷 月 最 低温 度 (bio6)〉 的 范围 是 -6.17~3 ‘C。6 个 主导 变量 的 
变异 系数 为 17.20% 一 79.37%, 最 冷 月 份 最 低温 (bio6 ) 的 变异 系数 最 大 , 最 湿 季 降雨 量 (bio16) 
的 变异 系数 最 小 。 综合 各 气候 变量 特点 来 看 , 石山 曹 苦 属 适宜 分 布 区 整体 上 具有 湿热 的 气候 
特点 GES) 。 
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HA 


降雨 量 Precipitation of the wettest quarter (mm) 最 热 季 度 平均 降雨 量 Mean precipitation of the warmest quarter (mm) 


A. 温度 变化 方差 (bio4) ; B. 海拔 CalO ; C. 最 冷 月 最 低温 度 (bio6) ; D. 最 干 月 降雨 量 (biol4) ; 
E. 最 湿 季 降雨 量 (biol6) ; F. 最 热 季 度 平均 降雨 量 (bio18) 。 
A. SD of temperature seasonality ( bio4); B. Altitude (alt); C. Min temperature of the coldest month (bi06); D. 


Precipitation of the driest month (bio14); E. Precipitation of the wettest quarter (bio16); F. Mean precipitation of 


the warmest quarter (bio18). 


DS 


4 环境 变量 的 响应 曲线 和 物种 存在 概率 
Fig. 4 Response curve of environmental variables and probability of Petrocodon species 
distribution 


K 5 Ali EE a Eh Ai Kaa E IR Sti 
Table 5 Descriptive statistics of environmental variables in the distribution regions of 


Petrocodon 
K o. E 95% E fi EX [H] 
环境 变量 均值 。” 标准 差 ”变异 系数 m 最 小 值 RK 
95% confidence interval 
下 限 ER 
Environmental Standard 
l Mean u CV (95) Minimum Maximum 
variable deviation Lower bound Upper bound 
alt 710.42 413.71 58.23 649.40 771.45 70 2084 
biol4 28.26 16.07 56.88 25.89 30.63 1 54 
biol6 741.52 127.51 17.20 722.71 760.33 406 1 207 
biol8 653.90 164.14 25.10 629.69 678.11 185 1 207 


bio6 -1.47 4.31 79.37 70.144. —0.087 76.17 3 


95% 4 fa [X [8] 


环境 变量 均值 PREZ 变异 系数 最 小 值 最 大 值 
95% confidence interval 
下 限 ER 
Environmental Standard 
. Mean 2 CV (95) Minimum Maximum 
variable deviation Lower bound Upper bound 
bio4 596 156 26.22 573 619 98 920 


2.3 Aue ER 2A ae EE AREE 
对 120 个 地 理 信 息 在 ArcGis PHETEM, ARER, GUA BW es dz EXE 
续 分 布 在 中 国 西 南部 至 中 南 半 岛 北部 ， 零 散 分 布 于 华南 沿海 和 华中 地 区 。 预 测 分 布 图 5. 图 
6 表明 ，120 个 石山 昔 苦 属 分 布 记录 点 中 有 103 个 点 被 包含 在 潜在 的 适宜 区 内 ， 表 明 预 测 结 
果 与 实际 分 布 基本 保持 一 致 。 在 当前 气候 情景 下 ， 总 适宜 栖息 地 面积 为 1.87x107 km?， 占 中 
国 和 中 南 半 岛 陆地 面积 的 11.66%; 高 适宜 栖息 地 面积 为 9.980x105 km?， 占 研究 区 总 陆地 面 
一 积 的 6.21%; 最 适宜 栖息 地 面积 为 8.76x105 km?， 占 总 面积 的 5.45%。 结 果 表 明 ， 当 前 石山 
iJ 营 苔 属 总 适宜 生境 范围 为 104° 一 112* E, 18°—32°N, JUPE P gi E ERE CE, 
零散 分 布 华南 沿海 和 华中 地 区 ， 块 状 分 布 于 缅甸 北部 。 
2.4 气候 变化 下 石山 曹 营 属 分 布 格局 和 面积 的 时 空 变化 
根据 当前 的 地 理 分 布 记录 和 环境 数据 ， 绘 制 了 石山 曹 苔 属 潜 在 分 布 区 预测 图 (图 5， 图 
-— 6) 。 末 次 间 冰 期 适宜 生境 面积 占 总 面积 的 10.73%， 在 当前 时 期 总 适 生 区 面积 的 基础 上 新 增 
v FEN 2.26%, EREN 3.9690, BPD TEC. SHIT HIG TAR EL 11.66%, # 
S 占 比 高 于 末次 间 冰 期 但 当前 适宜 生境 的 分 布 中 心 相 比 末次 间 冰 期 向 高 纬度 山区 扩散 ,并且 
c 适宜 区 破碎 化 增强 。 比 较 末 次 冰 盛 期 与 当前 潜在 分 布 区 的 格局 表明 , 末次 冰 盛 期 几乎 没有 该 
e 属 植物 适宜 的 分 布 区 ， 形 失 率 高 达 11.71%， 远 大 于 新 增 率 ， 表 明 气 候 寒冷 干燥 不 利于 石山 
营 苔 属 的 生长 。 全 新 世 中 期 与 当前 相 比 ,该 属 植物 的 适宜 生境 轮 廊 与 当前 相似 , 适 生 区 面积 
占 比 为 10.33%， 新 增 率 0.86%, TERA 3.26%， 新 增 率 小 于 丧失 率 ， 表 明 该 时 期 气候 条 件 仍 
(U RAAT AWE BBE. 在 历史 气候 条 件 下 ， 自 末次 间 冰 期 以 来 , 云贵 高 原 至 越南 北部 
= 高 度 适 生 区 均 明显 缩减 与 破碎 化 。 
= 在 未 来 气候 情景 下 (图 7 ) ， 于 2050 FALE KERERE LL (14.66%) 达到 峰 
值 ， 新 增 率 2.6196, ERK 0.04%， 新 增 区 域 位 于 当前 时 期 适宜 分 布 区 中 心 的 东部 和 北部 。 
2070 年 适宜 生境 占 比 为 14.17%， 新 增 率 2.7296, WRZ 0.24%。 未 来 两 个 时 期 新 增 区 域 均 
位 于 现代 时 期 适宜 生境 的 北部 和 东部 , 这 是 物种 对 未 来 气候 变化 反应 的 敏感 区 域 , 这 些 区 域 
应 为 重要 保护 区 域 。 同时 ,表明 在 2050 气候 变 暖 的 环境 条 件 下 ， 高 度 适 生 区 呈现 集中 分 布 ， 
并 且 石 山 营 苔 属 植 物 的 分 布 在 一 定 的 温度 升 高 范围 内 处 于 积极 状态 ， 直 到 2070 将 会 随 着 温 
度 的 升 高 适宜 区 面积 有 所 降低 。 
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A. 末次 间 冰 期 (LIG, IES% 120 000—140 000 年 ) ; B. 末次 冰 盛 期 (LGM, K4 22000 年 前 ) ; C. 全 
新 世 中 期 C(MH， 大 约 6 000 年 前 ) ; D. 当前 (1950 一 2000 年 ) ; E.2050 年 (2041—2060 年 ) ; F. 2070 


A. Last interglacial (LIG, about 120 000-140 000 years BP); B. Last glacial maximum (LGM, about 22 000 years 
ago); C. Mid-holocene (MH, about 6 000 years ago); D. Current (1950—2000); E. 2050 (average for 2041—2060); 


F. 2070 (average for 2061—2080). 
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Fig. 5 Potential geographic distribution of Petrocodon 
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K 6 不 同时 期 石山 营 苔 属 适 生 区 分 布 面积 的 百分比 (单位 : %) 
Table 6 Percentage of suitable area of Petrocodon at different periods (Unit: %) 
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时 期 不 适 生 区 低 适 生 区 高 适 生 区 最 适 生 区 总 适 生 区 

Period Unsuitable Lowly Highly Most suitable All suitable 
area suitable area suitable area area area 
末次 间 冰 期 Last interglacial (LIG) 81.55 7.52 4.08 6.65 10.73 
末次 冰 盛 期 ” Last glacial maximum (LGM ) 83.62 13.40 2.56 0.42 2.99 
全 新 世 中 期 Mid-holocene (MH) 80.02 9.66 5.96 4.37 10.33 
当前 Current 80.27 8.07 6.21 5.45 11.66 
2050 年 2050 78.47 7.24 5.19 9.09 14.29 
2070 年 2070 78.74 7.09 5.25 8.92 14.17 
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A. 当前 相对 于 LIG; B. 当前 相对 于 末次 冰 盛 期 ，C， 当 前 相对 于 全 新 世 中 期 ，D， 当前 相对 于 2050 4E; E. 
当前 相对 于 2070 4E; F. 石山 蕊 苔 属 分 布 记录 。 不 同时 期 缩写 含义 详 见 表 6。 
A. Current relative to LIG; B. Current relative to LGM; C. Current relative to MH; D. Current relative to 2050; E. 


Current relative to 2070; F. Distribution records of Petrocodon. See Table 6 for abbreviations of different periods. 
图 6 DS RISE UTER a LL E S Dg EE eA, 


Fig. 6 Changes in potential suitable habitats for Petrocodon in historical and future climatic 


scenarios 
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不 同时 期 缩写 详 见 表 6。 下 同 
See Table 6 for abbreviations of different periods. The same below. 
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Fig. 7 Change rate ( % ) of suitable areas of Petrocodon under different scenarios 


2.5 气候 变化 背景 下 石山 曹 若 属 潜在 分 布 中 心 的 迁移 路 线 

末次 间 冰 期 适宜 区 质心 (图 8) 位 于 广西 永福 县 (110.10° E, 24.69? ND ， 末 次 间 冰 期 
向 西北 迁移 到 贵州 丹 寨 县 (107.94° E, 26.16°N) ， 末 次 冰 盛 期 向 东南 部 迁移 至 广西 三 江 县 
南部 (109.60° E, 25.54? N)， 全 新 世 中 期 再 向 东北 方向 迁移 至 湖南 省 通道 县 东部 (109.95° 
E, 26.24° N) ， 之 后 向 东部 迁移 至 湖南 城 步 县 西北 部 (110.29* E, 26.05? ND , 2070 年 
质心 点 向 东南 迁移 湖南 城 步 县 南部 (110.19° E, 26.30°N) 。 末 次 间 冰 期 相 比 当前 向 南 迁移 
179 km; 末次 冰 盛 期 相 比 当前 向 西 迁移 194 km; 全 新 世 中 期 相 比 当前 向 西南 移动 87 kmo 
此 外 ，2050 年 向 东北 迁移 25 km、2070 年 东南 方向 迁移 40 km。 自 末次 间 冰 期 以 来 ， 石 山 
营 苔 属 适 宜生 境 总 体 向 高 纬度 地 区 迁移 , 迁移 至 湖南 省 南部 为 最 北 界 , 并 且 各 时 期 迁移 距离 
存在 较 大 差异 。 
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A. 研究 区 域 ，B. 迁移 路 线 。 
A. The study area; B. Migration route. 


图 8 不 同 气候 情景 下 石山 苞谷 属 潜在 分 布 中 心 的 迁移 路 径 


Fig. 8 Migration paths of potential distribution centers of Petrocodon under different climatic 


scenarios 
3 讨论 与 结论 
3.1 模型 精度 分 析 
MaxEnt 模型 的 参数 默认 设置 是 通过 对 来 自 6 个 地 区 的 226 种 物种 的 不 同 数 据 集 进行 测 


试 而 生成 的 (Phillips & Dudik, 2008) ， 应 用 ROC 曲线 评价 模型 精度 。 但 是 对 于 特定 物种 的 
研究 ， 使 用 专门 的 模型 参数 设置 比 使 用 默认 设置 更 有 效 CAnderson & Gonzalez, 2011) ， 可 
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以 有 效 的 避免 过 度 拟 合 (Zhu & Qiao, 2016) 。 优 化 的 MaxEnt 模型 预测 结果 得 到 平滑 的 响应 
曲线 且 AUC 宇 0.96， 表 明 微 调 模 型 可 以 合理 地 反映 该 属 植物 地 理 分 布 对 气候 变量 变化 的 响 
应 〈 朱 耿 平 和 乔 慧 捷 ，2016) ， 并 可 以 准确 预测 其 潜在 分 布 。 表 明 此 模型 用 于 石山 昔 苔 属 潜 
在 分 布 区 的 预测 可 信 度 高 。 
3.2 石山 童 昔 属 潜在 适 生 区 现状 分 析 
MaxEnt 模型 预测 结果 表示 石山 营 苔 属 分 布 的 最 大 可 能 性 ， 不 能 准确 得 出 物种 现实 的 分 
布 范围 ， 而 是 物种 的 潜在 分 布 区 。 我 们 的 结果 表明 : 从 末次 间 冰 期 (LIG) 开始 至 2070 年 
经 过 6 个 时 期 , 石山 营 苔 属 潜在 区 从 中 国 西南 部 连续 至 越南 北部 , 零散 分 布 中国 南 部 和 中 部 
及 块 状 分 布 缅甸 北部 。 因 此, 中 国 西南 部 至 越南 北部 的 地 形 地 和 瑶 优 势 极 有 可 能 是 石山 莒 苦 属 
避难 的 主要 场所 。 这 些 地 区 为 石山 鞋 苦 属 适 生 区 ,也 是 石山 莒 苦 属 种 质 资源 热点 地 区 ， 是 未 
来 石山 曹 苦 属 重点 保护 的 区 域 。 
如 前 文 所 述 , 季风 是 石山 莒 若 属 从 典型 热带 地 区 扩散 至 热带 北 缘 的 中 国 西南 地 区 并 存活 
下 来 的 先决 条 件 (Ren, 2015) ， 且 高 山河 谷 与 喀斯特 地 貌 是 该 属 植物 的 主要 分 布 区 域 。 
究 区 是 亚洲 热带 的 北 缘 以 及 南亚 季风 、 东 亚 季 风 和 西北 太平 洋 季 风 的 又 加 影响 区 , 也 是 高 
河谷 与 喀斯特 地 貌 集中 分 布 区 〈 姜 超 等 ，2017) ,该 区 域 局 域 高 度 异 质 性 的 地 形 地 貌 造 成 多 
> 样 的 生境 有 利于 石山 曹 苦 属 的 生长 与 繁殖 , ET A EA LU E A IE INE SES 
N 性 干 湿 气 候 、 石 灰 岩 高 钙 土 壤 等 特殊 环境 。 从 季节 角度 分 析 发 现 ， 南 亚 季风 在 5 一 7 月 强度 
最 大 《〈 周 浙 昆 等 ，2017) ， 东 亚 季 风 6 一 8 月 强度 最 大 《〈 吴 国雄 等 ，2013) ， 西 北 太 平 洋 季 
风 5 一 9 月 强度 最 大 〈 黎 舌 等 ，2014) ， 在 季风 作用 最 强 的 阶段 正好 为 该 属 植物 生长 发 育 季 
之 节 。 因 此 ， 在 石山 莒 昔 属 的 引种 驯化 实践 中 ， 在 选择 引种 植 基地 时 着 重 考虑 当地 季风 情况 。 
一 3.3 环境 变量 与 石山 莒 苦 属 潜 在 地 理 分 布 变化 的 关系 
Lam 由 于 石山 营 苔 属 本 身 的 生理 特性 和 所 在 生境 的 特殊 性 , 使 其 仅 分 布 于 特定 地 区 或 某 种 局 
部 特殊 生境 。 在 不 同时 期 内 ,影响 石山 营 苔 属 分 布 的 土壤 变量 和 地 形 地 貌 基本 相同 , 不 会 有 
较 大 波动 , 而 致使 石山 昔 苔 属 分 布 区 变化 的 可 能 是 由 于 气候 变量 的 剧烈 波动 以 及 人 类 活动 的 
强烈 干扰 。 当 前 ， 制 约 石山 营 芭 属地 理 分 布 的 主导 环境 变量 为 最 干 月 降雨 量 (bio14) 、 最 
热 季 度 平均 降雨 量 (bio18) 、 最 湿 季 降雨 量 (bio16) 、 温 度 变化 方差 (bio4) 、 最 冷 月 最 
IEEE (bio6) 、 海 拔 (alt) 。 可 见 ， 温 度 和 降水 作为 重要 的 环境 变量 显著 影响 着 石山 莒 若 
属 的 生理 特征 、 分 布 范 围 、 多 样 性 和 丰富 度 等 ; 海拔 变量 也 具有 显著 影响 作用 是 通过 重新 分 
配 降水 和 温度 来 实现 (Lambert et al., 2010; Wang et al., 2014) 。 

ALESE 5 一 9 月 开花 、10 一 12 月 形成 成 熟 藉 果 ， 整 个 营养 生长 和 开花 结果 阶段 正好 
与 研究 区 最 暖 季度 重合 。 因 此 ， 干 旱 和 寒冷 季度 可 能 会 阻止 该 属 植物 的 繁衍 。 同 时 ， 温 度 变 
化 方差 (bio4) 和 最 冷 月 最 低温 (bio6) 为 主导 环境 变量 合理 性 在 于 一 定 温 差 和 最 低温 是 植 
物 的 生长 、 发 育 和 开花 最 重要 调控 因子 之 一 (Khodorova et al, 2014) ， 但 是 ， 低 于 -6.17 °C 
极端 低温 不 利于 石山 营 苔 属 的 生长 ， 从 而 影响 其 在 地 理 尺度 上 的 分 布 格局 。 因 为 石山 营 苔 属 
从 热带 或 湿热 地 区 进化 而 来 ， 可 能 不 具备 抵御 寒冷 的 机 制 (Li, 1996〉。 因 此 ， 在 石山 营 苔 
属 的 引种 栽培 中 , 在 选取 作为 选 育 优质 种 源 的 种 植 基地 时 要 注意 当地 温度 变化 方差 和 最 冷 月 
最 低温 ， 应 优先 选择 最 冷 月 最 低温 为 -6.17~3 'C， 温 度 变化 方差 是 98 一 920 的 区 域 。 

此 外 ,Liu 等 (2017) 研究 表明 降水 和 海拔 都 不 会 显著 影响 苦 曹 苦 科 植物 地 理 分 布 模式 ， 
但 是 本 研究 表明 ， 最 干 月 降雨 量 (bio14) 、 最 热 季 度 平均 降雨 量 (bio18) 、 最 湿 季 降雨 量 
(biol6) 和 海拔 (alt) 对 石山 曹 苦 属 的 地 理 分 布 有 着 显著 的 影响 ， 可 能 是 由 于 石山 莒 苦 属 
主要 生长 在 喀斯特 岩 壁 上 ， 最 热 季 度 平 均 降 雨量 (185 一 1 207 mm) 可 加 速 石灰 岩 的 矿 化 速 
率 ， 可 为 其 生长 发 育 提供 了 更 多 的 养分 ( 徐 翔 等 ，2018) 。 石山 营 苔 属 的 存在 概率 随 海拔 的 
上 升 而 减少 的 趋势 ,， 且 峰值 出 现在 中 低 海 拔 段 (0 一 1 500 m) ,与 实际 分 布 海拔 C120—1 700 
m) 相 比 预测 分 布 区 范围 更 广 ， 但 是 超过 1500 m 的 高 海拔 区 域 不 利于 大 多 数 石山 莒 营 属 的 
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生长 。 正如 水 分 能 量 假说 所 认为 , 水 分 和 热量 的 共同 作用 影响 着 植物 物种 分 布 的 格局 CA Ten, 
2002; Hawkins et al., 2003; Bradford et al., 2003; 王 志恒 等 ，2004a) 。 所 以 ， 我 们 推测 ， 高 海 
拨 区 域 存在 概率 低 主要 是 由 于 低温 环境 降低 了 石灰 岩 的 矿 化 速率 ; 喀斯特 地 貌 特 殊 的 双 层 地 
质 结构 ,致使 水 分 迅速 流失 ， 从 而 形成 了 干旱 低温 缺 水 的 喀斯特 地 貌 环 境 不 利于 石山 苔 苔 属 
的 生存 (Hawkins et al., 2003; 王 志 恒 等 ，2004b) 。 因 此 ， 在 石山 莒 苔 属 未 来 的 生产 实践 中 ， 
要 注意 当地 降水 量 ， 尤 其 要 重点 关注 最 热 季 度 平均 降雨 量 的 适宜 取 值 在 185—1 207 mm、 最 
干 月 降雨 量 适 宜 取 值 在 1—54 mm、 最 湿 季 降雨 量 适 宜 取 值 在 406~ 1 207 mm 的 地 域 。 

综 上 所 述 , 各 项 环境 变量 通过 控制 石山 曹 若 属 的 生长 和 发 育 规律 ， 从 而 进一步 影响 其 分 
布 格局 。 能 量变 量 、 海 拔 变 量 、 水 分 变量 在 塑造 石山 营 苔 属 生态 适应 方面 发 挥 着 主要 作用 ， 
j 十 壤 变 量 对 其 影响 较 小 , 可 能 是 因为 石山 营 昔 属 多 在 喀斯特 地 区 岩石 表面 半 附 生 至 附 生 生 
活 , 在 大 尺度 区 域 条 件 下 土壤 的 理化 性 质 不 足以 代表 石山 曹 苦 属 生长 发 育 的 土壤 机 制 。 其 次 ， 
在 进行 地 理 分 布 数据 统计 以 及 查阅 相关 研究 发 现 , 石山 曹 苔 属 植物 地 理 分 布 信息 均 来 自 喀 斯 
特地 貌 〈《 除 星 苯 石 山 童 苦 ) ， 因 此 在 石山 曹 苦 属 〈 除 星 苯 石 山 曹 苦 ， 该 种 目前 仅 有 已 知 其 分 
布 广西 资源 县 和 湖南 新 宁县 两 省 区 县 交界 的 丹 起 地 貌 上 )〉 所 有 种 均 在 相同 的 喀斯特 地 貌 中 。 
在 此 条 件 下 , 能 够 反映 喀斯特 地 貌 缘 景 的 土壤 变量 的 差异 不 能 从 该 属 植 物种 间 的 地 理 分 布 中 
~> 表现 出 来 。 

N 3.4 A l Ë E REE K 8 Ae AA AR LESER SZ 
根据 模型 结果 ， 从 末次 间 冰 期 开始 到 未 来 气候 (2070 年) 6 种 气候 情境 中 ,石山 营 

在 中 国 及 中 南 半 岛 的 适宜 生境 面积 在 末次 冰 盛 期 断崖 式 降低 、 全 新 世 中 期 成 倍增 加 、 之 
e 慢 增 加 ， 直 至 2050 年 停止 增加 、 于 2070 年 出 现 缓慢 降低 趋势 ， 且 整体 有 向 高 纬度 地 区 
r 的 趋势 。 根 据 气 候 稳定 性 假说 ， 气 候 波 动 会 导致 物种 濒危 ， 而 稳定 的 气候 往往 会 增加 一 个 地 
Al 区 的 物种 丰富 度 〈Stevens, 1989) 。 在 来 次 间 冰期 间 ， 东 亚 季风 增强 、 降 雨 增 加 ， 全 球 温暖 
程度 与 当今 的 增 暖 相当 ， 适 生 区 分 布 接近 当代 分 布 格局 。 但 是 适 生 区 面积 相 比 当前 较 少 ,可 
能 是 因为 在 末次 间 冰 期 植被 的 反馈 作用 对 华南 地 区 降水 的 影响 不 显著 ( 张 琼 和 陈 婕 ,2020)。 
在 末次 冰 盛 期 ,气候 寒冷 干燥 、 波 动 大 ,不 利于 石山 营 苔 属 的 生长 发 育 ， 导 致 在 该 时 期 适 和 9 
区 面积 丧失 率 高 达 11.71%。 在 末次 冰 盛 期 之 后 全 球 气候 进入 了 温度 急剧 上 升 的 冰 消 期 ， 如 
球 气 候 逐 步 进 入 了 相对 稳定 的 全 新 世 暖 期 〈 况 雪 源 等 ，2021) ， 该 时 期 适宜 区 的 北部 和 东北 
部 急速 向 外 缘 扩 增 ， 适 生 区 面积 占 比 达到 10.33%。 当 前 时 期 至 2050 年 ， 温 度 的 逐渐 升 高 促 
使 了 石山 童 苦 属 适 生 区 的 扩展 ， 使 其 适 生 区 面积 占 比 达到 最 大 值 14.29%， 说 明 一 定 范 围 的 
气候 变 暖 有 利于 石山 营 苔 属 的 扩散 。 但 是 ， 在 2070 年 ， 适 宜生 境 相对 有 所 降低 ， 说 明 全 球 
气候 变 暧 也 可 能 致使 该 属 物 种 实际 生态 位 变 罕 ,因为 石山 昔 苔 属 对 石灰 岩 土壤 基质 有 特殊 要 
5K (Liu etal., 2017) ， 它 们 难以 跟踪 适宜 气候 的 北 移 。 

总 体 看 来 , 石山 曹 苦 属 在 中 国 及 中 南 半 岛 相对 稳定 的 高 适 生 区 主要 为 中 国 西南 部 (广西 、 
贵州 、 云 南 东南 部 ) 和 越南 北部 。 在 这 些 地 区 加 强 保护 力度 ， 对 自然 保护 区 进行 “禁止 滥 采 
Wits” 法 律 宣 传 。 并 且 ， 石 山 昔 苔 属 适宜 分 布 区 整体 向 高 纬度 地 区 移动 的 现象 ， 这 与 前 人 
究 结果 保持 一 致 〈《Sekercioglu et al., 2008; 方 精 云 等 ，2018) 。 适 宜生 境 中 的 收缩 和 扩张 位 
于 适宜 区 域 的 北部 和 东北 部 ， 这 是 物种 对 气候 变化 反应 的 敏感 区 域 CThuiller et al., 2005; 
Diamond et al., 2011) ， 这 些 区 域 应 该 加 强 监 测 力度 ， 针 对 性 制定 更 好 的 指导 保护 策略 。 

£x E, 鉴于 石山 营 苔 属 大 多 数 面临 濒危 的 现状 和 具有 较 高 的 观赏 价值 , 建议 基于 石山 蔚 
苦 属 自然 群落 在 高 适 生 区 实施 系统 详细 调查 ， 摸 清 其 致 濒 机 制 ， 基于 潜在 适宜 分 布 区 ， 开 展 
石山 鞋 营 属 引种 、 开 发 利用 研究 ， 为 园林 观赏 植物 的 应 用 做 好 基础 准备 。 

BR 谢 在 这 里 非常 感谢 中 国 科 学 院 昆明 植物 研究 所 税 玉民 研究 员 提供 免 儿 风 叶 石山 
EE (Petrocodon ainsliifolius W.H. Chen & Y.M. Shui) 精 确 的 经 纬度 地 理 信息 ; 也 非常 感谢 江 
西 省 中 国 科学 院 庐山 植物 园 的 助理 研究 员 张 涛 为 本 文 研究 提供 中 国 行政 区 矢量 图 。 
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